Model-Independent Test of General Relativity: 
An Extended post-Einsteinian Framework with Complete Polarization Content 
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We develop a model-independent test of General Relativity that allows for the constraint of 
the gravitational wave (GW) polarization content with GW detections of binary compact object 
inspirals. We first consider three modified gravity theories (Brans-Dicke theory, Rosen's theory and 
Lightman-Lee theory) and calculate the response function of ground-based detectors to gravitational 
waves in the inspiral phase. This allows us to see how additional polarizations predicted in these 
theories modify the General Relativistic prediction of the response function. We then consider 
general power-law modifications to the Hamiltonian and radiation-reaction force and study how 
these modify the time-domain and Fourier response function when all polarizations are present. From 
these general arguments and specific modified gravity examples, we infer an improved parameterized 
post-Einsteinian template family with complete polarization content. This family enhances General 
Relativity templates through the inclusion of new theory parameters, reducing to the former when 
these parameters acquire certain values, and recovering modified gravity predictions for other values, 
including all polarizations. We conclude by discussing detection strategies to constrain these new, 
polarization theory parameters by constructing certain null channels through the combination of 
output from multiple detectors. 

PACS numbers: 04.80.Cc,04.80.Nn,04.30.-w,04.50.Kd 
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I. INTRODUCTION 

Gravitational waves (GWs) will soon be detected by 
ground-based detectors, such as Advanced LIGO 
and Advanced Virgo [j, HI. These waves will provide 
invaluable information about the gravitational interac- 
tion in the so-called strong-field regime, where the non- 
linearity and strong dynamics of the Einstein equations 
play an important role. Although General Relativity 
(GR) has passed all Solar System and binary pulsar tests 
with flying colors, the strong-field regime remains mostly 
unexplored Q- For example, observations have not 
yet been able to confirm GR's no-hair theorems or the 
non-linear part of the Einstein equations in the GW gen- 
eration. 

One of the primary targets of these terrestrial detectors 
are waves generated in the late inspiral of compact ob- 
jects. Detectors that can observe signals as low as 10 Hz 
will be able to follow the binary inspiral from astronomi- 
cally small separations down to merger. For example, for 
a binary neutron star inspiral, ground-based detectors 
should be able to detect inspirals from initial separations 
of ^ 600mNs ~ 2500 km, where ttins is the total mass 
of the binary (here and in what follows we use geomet- 
ric units G = c = 1). In this late inspiral regime, the 
non-linearities of the field equations and the strong-field 
nature of GR are essential, while one can still employ 
perturbation theory to model the associated GWs until 
right before the compact objects plunge into each other 
and merge. 

One of the predictions of GR that one would wish to 
test is that the GW metric perturbation only possesses 
two propagating degrees of freedom. In a general theory 
of gravity, there are up to 6 degrees of freedom allowed. 



In GR, however, due to the structure of the field equa- 
tions, only 2 of them are physical, with the remaining 4 
being gauge degrees of freedom. Although ground-based 
detectors are expected to observe GWs at low signal-to- 
noise ratio, the simultaneous detection of such waves by 
multiple detectors should allow us to constrain the ex- 
istence of additional polarization modes (2 scalar and 2 
vectorial). The main question, of course, is precisely how 
to carry out such tests and the degree to which they will 
allow us to rule out or confirm additional polarization 
modes. 

Many modified gravity theories exist where all 6 de- 
grees of freedom are physical and must be accounted for 
when computing the GW response function. Perhaps 
the most well-known examples are scalar-tensor theories, 
where the presence of a scalar field leads to the exis- 
tence of dipole radiation and a scalar (so-called breath- 
ing) polarization mode, in addition to the two standard 
transverse-traceless modes of GR Vector-tensor 
theories usually predict the existence of preferred direc- 
tions and the excitation of vector modes p, while tensor- 
vector-scalar theories, such as TeVeS (9l.]lO|. predict the 
existence of all 6 polarization modes and Einstein- Aether 
theories [TT| - [l3j predict the existence of 5 polarization 
modes . Moreover, bimetric and stratified theories also 
predict the existence of all 6 polarization modes @, [l^ - 

Recently, a model-independent framework to test 
GR was developed: the parameterized post-Einsteinian 
(ppE) framework [13, [13 ■ In this scheme, one enhances 
the GR GW template family through the addition of new 
theory parameters 6; when these ppE parameters acquire 
certain values one recovers GR, while for other values one 
recovers predictions from modified gravity theories. As 
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a first step, Yunes and Pretorius [l7[ proposed a sim- 
ple ppE model that treated the inspiral phase in a post- 
Newtonian (PN) approximation, where the orbital veloc- 
ity is much less than the speed of light w/c ^ 1 , and that 
neglected the direct excitation of additional polarization 
modes in the GW response function. Such modes were 
only indirectly accounted for through possible modifica- 
tions to the orbital binding energy and its rate of change, 
which were captured through ppE parameters. 

In this paper, we relax this assumption and improve 
the ppE framework to allow for the direct presence of ad- 
ditional polarization modes in the response function. In 
order to achieve this, we first look at three modified grav- 
ity theories that predict the existence of additional po- 
larization modes, namely Brans-Dickc theory Rosen's 
bimctric theory [13, [HI and Lightman-Lee theory [l^ ^ . 
For each of these theories, we calculate the GW met- 
ric perturbation, we extract from this the polarization 
modes, construct the GW response function and Fourier 
transform the latter in the stationary-phase approxima- 
tion (SPA) [HHll. 

We then use physical arguments to predict the general 
functional form of the time-domain response function, 
assuming all polarization states are present. We begin 
by noting that the existence of additional modes leads 
to new terms in the time-domain response function that 
are proportional to the t = (0, 1, 2) harmonics of the or- 
bital phase. When this is Fourier transformed, however, 
the ^ = harmonic docs not possess a stationary point, 
and thus, it is subdominant relative to the I = (1,2) 
modes. When Fourier transforming, we allow for a para- 
metric deformation of the Hamiltonian and the radiation- 
reaction force through the addition of relative power-law 
corrections. As Yunes and Pretorius found [17[, power- 
law modifications to the binding energy and energy flux 
introduce power-law modifications to the Fourier ampli- 
tude and phase. The inclusion of additional polarizations 
introduces a new term that is proportional to the £ = 1 
harmonic of the orbital phase and that was not included 
in the original ppE scheme. 

With these general arguments and solutions from mod- 
ified gravity theories at hand, we generalize the ppE 
framework to allow for the remaining four polarization 
modes in the response function. Given a single detector, 
we find that to parameterize only the £ = 2 harmonic one 
requires 4 ppE parameters for the most general case, in 
agreement with Yunes and Pretorius [TtI ; if one wishes to 
include all polarization modes, one needs a total of 9 ppE 
parameters instead. However, if one restricts attention 
only to power-law modifications to the Hamiltonian and 
radiation-reaction force in the generation of GWs, then 
for single-detectors one needs only 2 ppE parameters to 



^ We choose these theories as toy models because substantial work 
has already been done to analyze their properties [Tgl. |20|| . Once 
similar work is done for other theories, the analysis of this paper 
can be extended to include them. 



parameterize the £ = 2 harmonic and 4 ppE parameters 
to parameterize all polarizations. 

Part of our work is similar to a paper by Arun that re- 
cently appeared in the literature j24| . In that paper, he 
first defines the dipolar mode of the GW as the £ = 1 har- 
monic component of the wave. Then, assuming a general 
functional form for this component and a GR quadrupole 
plus dipolc frequency evolution, he calculates the Fourier 
transform of the response function. He finds that the 
Fourier response can be parametrized by 2 parameters: 
one that measures the amplitude of the dipolc mode rel- 
ative to the leading quadrupole mode; and one that mea- 
sures the relative strength of dipole emission. The reason 
that Arun manages to parameterize the waveform with 
a smaller number of parameter is that he neglects cor- 
rections to the conservative part of the dynamics and he 
restricts attention to dipole emission only. Our paper 
generalizes his results, allowing for generic deformations 
in the Hamiltonian and radiation-reaction force, while 
simultaneously accounting for all polarization modes. 

Using a single spherical [2^ or truncated icosohe- 
dral [2^ acoustic detector it is possible to extract and 
separate [27| the various polarization modes. In contrast, 
a standard two-arm interfcrometric detector only pro- 
vides a single projection of the polarization state, and it 
takes multiple independent projections from an array of 
detectors to fully constrain the polarization content. We 
construct a ppE model that could be used with multiple 
interferometers by restoring the beam-pattern dependen- 
cies of the ppE amplitude. This forces us to break-up the 
overall amplitude into several terms, each with new ppE 
theory parameters, thus enlarging the parameter space 
to 18 theory parameters in the general case and 10 for 
power-law corrections. With this at hand, we then dis- 
cuss a new multiple detector strategy to extract the ad- 
ditional polarization modes through the construction of 
null streams, ie. combinations that isolate specific polar- 
ization modes. One such null stream has the property 
that if GR were right, the stream would be free of gravi- 
tational wave energy for any designated sky location. If 
each detector output is thought of as a vector in signal 
space, then the construction of null channels reduces to 
the continuous projection of this signal in directions or- 
thogonal to the 2 GR polarization modes. In general, we 
find that at least 3 detectors are necessary to constrain 
additional modes, 4 detectors are necessary to constrain 
the vector modes and with 6 detectors it is possible to 
construct streams that arc null in every theory of gravity. 

The rest of this paper is organized as follows: Sec- 
tion [H] describes how the gravitational response func- 
tion is computed in a generic modified gravity theory 
and summarizes the SPA. Section IIIII discusses Brans- 
Dicke theory, Rosen's theory and Lightman-Lee's theory, 
computing the response function and its Fourier trans- 
form in the SPA for each of them. Section IIVI calculates 
the complete waveform for generic deviations in the sys- 
tems binding energy and balance law. Section |V] uses 
the results from the previous sections to construct a gen- 
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eralized ppE framework that allows for the existenee of 
additional polarizations modes. Section |Vl] discusses the 
construction of null streams for generic waveforms with 6 
polarization modes. Section [VIII concludes and discusses 
possible future research directions. 

We follow here the conventions of Misner, Thorne 
and Wheeler Greek letters stand for spacetime 

time indices, while Latin letters in the middle of the al- 
phabet fc,...) stand for spatial indices only; com- 
mas in index lists stand for partial derivatives and 
semi-colons for covariant derivatives; parentheses and 
square brackets in index lists stand for symmetrization 
and anti-symmetrization respectively such as ^(a^) = 
{l/2){Aap + AjBc) and {l/2){Aafs - Afsa); the 

metric is denoted via 17^,^ with signature (— , -f ); 
the Einstein summation convention and geometric units 
with G = c = 1 is assumed, unless otherwise specified. 



II. FOURIER TRANSFORM OF A GENERIC 
RESPONSE FUNCTION 

In this section, we construct the response function for 
a generic modified gravity theory, given the GW metric 
perturbation. We then briefly explain how the SPA to 
the Fourier transform of this response can be calculated. 
We conclude the section by showing how the algorithm 
works in the standard GR limit. 



A. Polarizations from the Metric Perturbation 

In this subsection, we mainly follow Will 0,0, [111. The 
response function of a detector to a wave with all possible 
polarizations is 

h{t) = F+h++F,,h''+F^^.h'"' + F^,M''+F^h'' + F^h'^ , (1) 

where F. are angular pattern functions and h' are wave- 
form polarizations. The former are given by [29| 



F, 



(1 + cos^ 9) cos 2ip cos 2(f) — cos 9 sin2V' sin 20 , 



(2) 



Fx = - ( 1 + cos 9) sin 2ip cos 20 -f- cos 9 cos 2?/^ sin 20 , 



F„ 

Fe 

F 



— sin 9{cos 9 cos 20 cos ip — sin 20 sin ■0) , 

— sin 6'(cos 9 cos 20 sin ip + sin 20 cos ip) , 

— - cos 20 sin^ 9 , 
■ cos 20 sin^ 9 . 
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The waveform polarizations can be computed from the 
contraction of certain basis vectors {efj,efj,ef,e^) (see 
e-g- [2^ 1131 noting that in [2^ (ef,ef), ie. the basis vec- 
tors orthogonal to N^, the unit vector pointing from the 



source to the detector, are denoted as 
waveform amplitudes A. , namely 



h'' = A,, 
h+ = etA: 



e^'-A'^ 



, 0i)) with the 

(8a) 
(8b) 
(8c) 



In these equations, A^, is the amplitude of the scalar 
breathing mode, is the amplitude of the scalar longitu- 
dinal mode. Ay are the amplitudes of the vectorial modes 
are the amplitudes of the transverse-traceless 

modes. 

The usual way to find the waveform amplitudes A. 
is to compute the linearized Riemann tensor evaluated 
with the trace-reversed metric perturbation [2^, [2^ . A 
more straightforward way to accomplish the same re- 
sult, however, is to construct operators that act on the 
trace-reversed metric perturbation directly and return 
the waveform amplitudes. In terms of these, the am- 
plitudes are given by 



A, 

A, 
At 



1 

2' 



(9a) 

(9b) 
(9c) 

(9d) 



where Pij = 6ij — NiNj is a projection operator orthog- 
onal to iV', a unit vector pointing from the source to the 
detector, while h^"^ is the trace-reversed metric pertur- 
bation and (5*^ is the Kronecker delta. 

One might wonder whether wc can reconstruct the full 
metric perturbation from the GW polarization modes in 
Eq. (jni). Notice, though, that Eq. ((9)) contains only 6 
degrees of freedom (1 in A^, 1 in 2 in A^, because it 
is transverse and 2 in A^t because it is transverse and 
traceless), while the full metric perturbation generically 
contains 10 degrees of freedom. Thus, for the inversion 
to be unique one must make a gauge choice, such as a 
pure traceless (yet not fully transverse) gauge. Doing so, 
the metric perturbation can be written as 



(3) 




(4) 




(5) 


w 


(6) 






(7) 





0, 

u 

3A 



A. 



:A, 



D \ 3 



2N(''A^ Al^ 



D 



D 



(10) 



Of course, such a metric reconstruction is unnecessary 
for our purposes because the observable is the response 
function and wc can project out the relevant degrees of 
freedom (those in Eq. ^) without making any gauge 
choice. 
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B. Stationary Phase Approximation 

In GW data analysis, one often works with the Fourier 
transform of the response function, which can be ob- 
tained analytically in the SPA j2l|-|23t . We here briefly 
review this method. The goal of the SPA is to compute 
the generalized Fourier integral 



(11) 



where the power n depends on how the frequency evolves. 
Such an integral vanishes for < n < 1 when the limits of 
integration arc ±oo, because, as the complex exponential 
oscillates, contributions from subsequent intervals can- 
cel out. Indeed, to first order we keep only the leading 
quadrupole emission for which n = 1/4. Therefore, we 
see that the i = harmonic of the orbital phase in the 
response function will contribute subdominantly to the 
SPA of the Fourier response. 



assuming that the response function h{t) is composed of 
a slowly varying amplitude A{t) and a rapidly varying 
phase £^{t), namely 

h{t) = A{t) (e*^*(*) + e-'^*(*)) , (12) 

with £ > the harmonic number and the orbital 
phase. 

Equation pT|) can be rewritten using Eq. (fT2)) as 



C. General Relativity Limit 

As an illustrative example, let us apply the above for- 
malism to the trace-reversed metric perturbation in GR. 
The trace-reversed metric perturbation for a two-body, 
quasi-circular orbit is given, to leading-order in the PN 
expansion, by (see e.g. [3l|) 



h{f) = / A{t) e2''*/*+*^*(*) + g27ri/t-j<?*(t) ^ j--[^3-j where the quadrupole moment is 

777 



The first term in square brackets does not have a station- 
ary point, ie. a value of t for which the derivative of the 
argument of the exponential vanishes, 2TTf + id^/dt = 0. 
Terms without a stationary point contribute subdomi- 
nantly to the generalized Fourier integral, and thus, they 
can be neglected by the Riemann-Lebesgue Lemma [2l| . 
The second term in Eq. does have a stationary 

point, which after Taylor expanding occurs when the 
derivative of the argument in the exponential vanishes, 
ic. F{to) f/e, where F{t) = $/(27r) is the orbital an- 
gular frequency. 

Wfth this at hand, the SPA of h{f) is [H 



(17) 



(18) 



Aito) 



(14) 



where the phase GW ^' is given by 

nF{t.)] = 2.j [£j^-iyF' + -. (15) 

In the rest of this paper, we will use these expressions to 
find an analytic representation of the Fourier transform 
of the response function in different modified gravity the- 
ories. 

As we will see in Sec. IIIIl the time-domain response 
function in a generic modified gravity theory will con- 
tain terms proportional to the harmonic of the orbital 
phase, as shown in Eq. (jl2[) . Terms in the Fourier trans- 
form of the response function proportional to the £ = 
harmonic are of the form 



fi is the reduced mass, m is the total mass, D is the 
distance to the source, r is the orbital separation and 
{x^,v^) are orbital trajectory and orbital velocity unit 
vectors. 

In order to explicitly calculate the waveform ampli- 
tudes via Eq. © we must first express (x%f)*) in the 
source system fc') 



a;' = cos <I> i* -I- sin $ f , 
= — sin $ cos $ f 



(19) 
(20) 



Choosing the coordinate system so that the vector from 
the source to the observer is on the y-z plane the vectors 
(iV*, ef,e\) and the polarization tensors e^, e^- are 



iV* = sin i f ^- cos L fc* 



cos L f — sin L fc* 
i(e-ef-Hefe-), 



(21) 
(22) 
(23) 

(24) 
(25) 



where we recall that $ is the orbital phase and l is the 
inclination angle. Applying the operators in Eq. we 
obtain 



(26) 



-dt. 



(16) 



where TT[-] is the transverse-traceless projection op- 
erator, and all other amplitudes are zero. We have 



5 



here imposed the Lorenz gauge condition /i'^^.^ = 0, 
which can be rewritten as /i^°.o = Njh^\o and /i™.o = 

Once we have the ampHtudes, we can compute the 
waveform polarization modes from Eq. ([8]) to obtain 

/i+ = -^cos2$(l + cos2 0, (27) 
Dr 

/i'^ = -^sin2$cost, (28) 
Dr 

and all other modes vanish. From Eq. ([T]), the response 
is then simply 

h'^'^it) = AoK^(27rA^F)2/3e-'2* + c.c. , (29) 



where 



A, 



(1 + cos^ b) — 2iFx cos l , 



(30) 



2$ is the GW phase and J\A = if/^ra is the chirp mass, 
with rj = mim2/m^ the symmetric mass ratio. We have 
here used Kepler's third law to simplify the final result, 
neglecting subdominant terms in the PN approximation. 

The SPA of the Fourier transform of this response func- 
tion can be computed following Sec. IIIBl Using the bal- 
ance law to relate the rate of change of binding energy 
to the GW luminosity, we can calculate the frequency 
evolution, which to leading-order in the PN expansion is 
given by 



dF _ 48 
'dt " 5ttM^ 



(2^AlF)"/3 [l + 0{u^)] . (31) 



With this in hand, we can now compute the well-known 
(restricted^ ie. leading-order in the amplitude) Fourier 
transform of the response function in the SPA, namely 



1/2 



A, 



OR — (^A^ ./)-'/ V-'*oH, (32) 



D 



where we have also introduced for future convenience 



4 256U^ ^ ^ \ n ^ n 

^ n— 



(33) 

although only the £ ~ 2 harmonic enters the GR wave- 
form. Here, (c™, Z™) are known PN coefficients that can 
be read for example from Eq. (3.18) in [33|, and we have 
defined the reduced ^-harmonic frequency 



Uf = 



/ 27rX/ 



1/3 



(34) 



such that U2 = (ttA^/)^/^, with / the GW frequency. 

Up until now, we have concentrate on the restricted 
PN approximation, but later on it will be important to 
determine whether modified gravity corrections to the 



Fourier amplitude are degenerate with PN amplitude cor- 
rections in GR. Amplitude corrections arise because the 
PN waveform contains an infinite number of higher har- 
monics, as one can see e.g. in Eq. (238) of [3l|. Therefore, 
the Fourier transform of such a waveform leads also to a 
sum of £ harmonic terms. The dominant one is the £ = 2 
mode, which was already described above in Eq. (|32l) . 
The next order terms are the £ = 1 and £ = 3 harmonics, 
which scale as [3J 



, 96/ """" D 
1/2 ^.2 



(35) 



(36) 



(1 3) 

where A^^ are amplitude factors that depend on dif- 
ferent combinations of the inclination and polarization 
angles. The key point here is that these terms enter at 
IPN order higher in the amplitude relative to the domi- 
nant £ = 2 mode, as can be established by looking at its 
j-5/6 frequency dependence. 



III. MODIFIED GRAVITY THEORIES 

A. Brans-Dicke Theory 

Brans Dickc theory Q is defined by the gravitational 
action (in Jordan frame) 



S 

"° " 16^ 



d'^x. 



-gm - -^cj.^'-cl,^^ - ^-^V] , (37) 



where g is the determinant of the metric g^j^, R is the 
Ricci scalar, (/> is a dynamical scalar field, y is a poten- 
tial for the scalar field and Wbd is a coupling constant. 
Usually, one sets the potential to zero, unless one is con- 
sidering massive Brans-Dicke theory [SJ] . Such a theory 
is a subset of scalar-tensor theories, where the coupling 
constant w(0) = cjbd = const. Variation of this action 
with respect to the metric and the scalar field leads to the 
modified field equations of the theory. Linearizing these 
field equations about a flat background ri^y , one can ob- 
tain evolution equations for hf^v We refer the interested 
reader to [10, H^l for further details. 

In Brans-Dicke theory, it is convenient to define the 
trace-reversed metric perturbation in terms of two other 
fields, a covariantly conserved tensor 9^'^ and the scalar 
field 6, to obtain 



(38) 



where (f>o is the asymptotic value of the scalar field at 
spatial infinity. The linearized field equations prescribe 
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the evolution of both cj) and 9^'^ ^ whose solution in the 
PN approximation is 

0'' = ^U-\dl^{V'v'~S:'x^), (39) 



— = -^S 
d^n D ' 



(40) 



where we have defined 



[{N ■ vY ~ {N ■ xY] 



/Gm 



-(Gr + 2A) 25* , 

r \ r 



1/2 



N ■ V 



(41) 



In these equations, we have also defined ^ = (2 + 
Wbd)"^ ^ (^BO for Wbd > 1, G = 1 - £,{si +S2- 2S1S2), 
5 = si — S2, where sa is the sensitivity of the Ath ob- 
ject, as defined in Brans Dicke theory, F = f — 2{miS2 + 
m2Si)/m and A = f — si — S2 + C'(C)- Clearly, Brans- 
Dicke theory reduces to GR in the Wbd — > 00 (or ^ — > 0) 
limit. 

With this at hand, we can now compute the polar- 
ization modes as in Sec. Ill Al Using the Lorenz gauge 
condition from [T^ = 0^, the waveform amplitudes 

are 



(42) 



where the longitudinal and vectorial modes vanish. The 
polarization modes are then 



4nS 




D 





















2Gfim 
Dr 



(43) 

cos2$(H-cos2i) > (44) 



Dr 



■ sin 2$ cos i , 



(45) 



where again the longitudinal and vectorial modes vanish. 

Putting all pieces together and simplifying expressions 
through the leading PN order expression for Kepler's 
third law 



2ttF 



Gt 



1/2 



(46) 



we find 



/l^°(t) = ABn^(2^7WF)2/3e-*2$ 



SBor/i/''^^(27rXF)i/3e-* 



CB^^{2nMFY/' 



■ c.c. , 



(47) 



2 The evolution equation for is □9'"' = -167rT'"', where t''" 
is a complicated extension of the Landau-Lifshitz pseudo-tensor. 
One can easily verify that this differential equation preserves the 
Lorenz gauge condition. 



where c.c. stands for the complex conjugate and where 
we have kept terms only linear in f . In these equations, 
we have also defined 

Abd = ^-F+(l + cos^ t) — 2ii^x COSi 



fcBD^+(l + cos^ i) + 2i ksoFx COS l + F-^ — sin^ l 



= ^GR + ^^BD.l ^ 

Sbd = e(-J^b5'sini) =eSBo,i 



Cbd — C 



-Y(r + 2A) 



(48) 
(49) 

(50) 



and fcBD = (1/2) + (2/3)(si -\- S2 - 2siS2). 

The Fourier integral of Eq. (|T7|) can be easily calcu- 
lated with the SPA, but this requires use of the orbital 
frequency evolution. The rate of change of the binding 
energy is determined by both dipole and quadrupole ra- 
diation, and to first order in ^ and to leading-order in the 
PN approximation, it is calculated in [l^ 



dE _ 8 ^i^TO^ 
H ~ "15 r4 



12G^ 



I 

12 



-G\S^ 
(51) 

When this is combined with the binding energy F ~ 
—Gm^/r and Eq. (j46p . one obtains the orbital frequency 
evolution 



dF_ _ 48 
'dt " 57r7W2 



{2nMF) 



11/3 



e {2TrMFy 

e (27rX^^)"/3 ( —r' - ^Bo 



(52) 



where recall that G 7^ 1 as given below Eq. (|^T|) . 

The total Fourier-transformed waveform in Brans- 
Dickc theory is then simply 

hnM) = hi'Mf) + hi^Uf), (53) 

where the Fourier transform of the first term in Eq. (j47p 
is 



;(2) 



96 ~d" 



A, 



1-e 



(- 

V24 



"2~ 



(54) 



and that of the second term is 

(55) 

We have here assumed that the second and third terms 
in Eq. (I52p arc much smaller than the first term. 
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ie. that Brans-Dickc theories introduces a small defor- 
mation away from GR. We have also here defined the 
^"^-harmonic Brans-Dicke Fourier phase 



where 



7168 



(56) 



with '^'il given in Eq. ^ (see also Eq. (Im)) ) and ut 
given by Eq. . Notice that the second term in Eq. ([M)) 
is of — IPN order relative to the first one, a typical signa- 
ture of dipole radiation. Such an amplitude correction is 
usually neglected, because GW interferometers are much 
more sensitive to the phase evolution. Of course, not all 
corrections to the amplitude and the phase will be mea- 
surable, and some of them can be degenerate with other 
system parameters, such as the luminosity distance or 
the chirp mass. Such issues will be discussed further in 
Sec. El 



B. Rosen's Theory 

Rosen's is an example of a bimetric theory [13, [Hi '■ a 
theory with a dynamical tensor gravitational field and a 
flat, non-dynamical metric, or prior geometry. Rosen's 
theory is defined by the gravitational action [^, [13, [lH 



5h 



1 
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where is the determinant of the flat, non-dynamical 
metric "q^u and the operator stands for a covariant 
derivative with respect to rj^^. Although the fleld equa- 
tions in this theory are quite different from Einstein's, 
it has a standard parametrized post-Newtonian (ppN) 
limit, with only the a2 ppN parameter different from 
its GR value. This parameter, however, has been con- 
strained to be less than 4 x 10~^ through observations of 
solar alignment with the ecliptic plane [^ [3] ■ 

Rosen's theory is of class IIq in the E{2) classiflca- 
tion [B, thus, not all six polarization modes are 

observer-independent. In theories of this class, all ob- 
servers agree on the magnitude of the longitudinal mode, 
but they disagree on the presence or absence of all other 
modes. This, however, does not mean that the other 
modes are not real or that they do not carry energy. It 
only implies that a spin-decomposition of the GWs is 
not invariant. Such frame dependence is irrelevant for 
our purposes, since the detector will be in a given frame, 
and thus, it will measure a certain number of polarization 
modes. See [1, [3] for more details on this theory. 

The lack of definite helicity in Rosen's theory is con- 
nected to the lack of positive definitencss in the sign of 
the emitted energy [36[- That is, for certain systems, 
Rosen's theory predicts dipolar radiation that pumps en- 
ergy into the system, leading to a total energy flux that 
is positive, instead of negative as in GR. In turn, this 
for example leads to an increase in the orbital period of 



binary pulsars with time |36l438| . Since binary pulsar 
observations are consistent with GR, Rosen's theory is 
today less appealing than in the 1970s. We here study it, 
not because we think of it as a particularly good candi- 
date to replace Einstein's theory, but as a toy model to 
determine how the GW response function is modifled in 
theories that allow for the existence of all gravitational 
polarization modes j^, [3] . 

The variation of the action and the linearization of the 
resulting field equations give the evolution of the trace- 
reversed metric perturbation [l^[33l, which can be solved 
to find 



l-(^.i)2]+ - 



^Lj^(7V.f))_a;J^(iV-£ 



— V V 



(58) 

(59) 
(60) 



where G is the difference in the self-gravitational binding 
energy per unit mass of the binary components: G = 
si/mi — 82/^2. Due to the specific characteristics of the 
theory, one cannot find a particular gauge to simplify the 
above equations. 

With this in hand, we can now compute the polariza- 
tions modes and the response function. Following the 
steps laid out in Sec. Ill Al we find 
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Using the modified Kepler's law 
2ttF 



m\i/2 



(62) 



where k„ = 1 — 4siS2/3, the time-domain response func- 
tion is 
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where we recall that c.c. stands for complex conjugate 
and where (A^, B^, C^) are functions of the angles: 



An = 



1 + COS^ i 

i cos i. — Fb 



sin^ L 



-F, sin L — F„A sin l ~ F,, 



sin2t 



-1/3 
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(66) 

Notice that in the G ^ limit, one does not recover GR, 
as described for example in [l^ [s^] ■ 

Before we can calculate the Fourier transform of the re- 
sponse function in the SPA for Rosen's theory, we must 
first calculate the orbital frequency evolution. The en- 
ergy evolution of the binary orbit due to GW emission 
is (to lea ding -order in the PN approximation and in G) 



given in [19 1 



dE 84^2^2 ^ 



dt 



15 r4 



20 fi'^m^G^ 



(67) 



and this leads to the orbital frequency evolution 



(68) 

where the binding energy is not modified. 

We can now calculate the Fourier transform in the 
SPA: 

hM) = hi'\f) + hi'\f), (69) 



where the transform of the first term in Eq. (|63p is 



84 D 

and that of the second term is 



,,/^:^(,A^/)-V6e-^*k^ (70) 



(71) 

We have here assumed that G ^ ^ and kept terms only 
to leading-order both in G and in the PN approximation. 
We have also here defined the ^-harmonic Rosen Fourier 
phase 
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4 ' " " ■ 224^5 

25^ fc-'/'g2^2/5 
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(72) 



where ui is given in Eq. p4p . As before, some modifica- 
tions are degenerate with system parameters, as we will 
see in Sec. IVl 



C. Lightman-Lee Theory 

Lightman-Lee theory [l^ is a bimetric theory of grav- 
ity, similar to Rosen's. This theory is controlled by the 
metric 5^1/, a dynamical gravitational tensor B^j^^ that is 
connected to and a flat, non-dynamical background 
metric @. The theory is defined by the gravitational 
action 

(73) 

where 77 is the trace of the background metric 77^^ (not to 
be confused with the symmetric mass ratio 77 introduced 
earlier). The spacetime metric is connected to the tensor 
Bf^„ via 



1 



1 

16 



(74) 
(75) 



which can be expanded for weak gravitational fields as 
S/xi- = V^Lv + with h^y = B^i, - B-q^^^/i. We refer the 
reader to @, 0] for more details on this theory. 

This theory is part of a wider class (sometimes called 
BSLL) of theories that are semi-conservative. In the par- 
ticular case of Lightman-Lee theory, there are no pre- 
ferred frame-effects and the ppN parameters reduce iden- 
tically to those of GR 0. This theory, however, is also 
of class IIq in the E{2) classification, just like Rosen's 
theory. As before, this implies the lack of definite helic- 
ity polarization states and the prediction of the existence 
of all polarization modes [36| . As in Rosen's theory, lack 
of helicity seems to lead to a lack of positive definiteness 
in the energy flux, which is in contradiction with binary 
pulsars [3^, Is^- Nonetheless, as in the case of Rosen's 
theory, we here take Lightman-Lee's as a toy model that 
allows us to see how the response function is modified in 
theories that predict all possible GW polarizations. 

The trace reversed metric perturbation is given in 



terms of by /i^^ 
in UM, we find that 



B„ 



lBrifj,ij. Using B as given 
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(77) 



h'^ = I± i^{4v'v^ - [2v^ - 2{N ■v)^ + 2 + 2{N ■ x)^]6'^} 
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(78) 



As in the case of Rosen's theory, one cannot choose 
a particular gauge to simplify expressions. We can then 
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extract the polarization modes as explained in Sec. Ill Al 
to obtain 
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With these modes, we can construct the response func- 
tion to find 
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where 
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In the derivation of the above equations we have used 
Kepler's law for this theory, which is identical to that of 
GR. 

Before wc can compute the Fourier transform of this 
response function in the SPA, we must first find the or- 
bital frequency evolution equation. The energy evolution 
of the binary's orbit is prescribed (to leading-order in the 
PN approximation and to leading-order in Q) by [l^ 



dE _ 84 /i^m^ 125 fi'^m'^g'^ 
lit ~ 15 rS ~9 



(83) 



Then, the orbital frequency evolution equation is simply 



dt 
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(84) 

The SPA of the Fourier transform of the response func- 
tion is 



(85) 



where the Fourier transform of the first term in Eq. (|79p 
is 



A^^i\l—'—(7rMf) ''"e 



(86) 



and that of the second term is 



D 

As before, we have here kept terms to leading-order in 
the PN approximation and in Q. Wc have also defined 
the ^^'^-harmonic Lightman-Lee Fourier phase 



^c- 271 ft. 
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where we recall that ui is given in Eq. (p4)) . As in the case 
of Rosen's theory, note that Lightman-Lee theory does 
not have a well-defined GR limit, ie. h^{f) 7^ ^gr(/) as 
— > 0. Also note that some modifications are degenerate 
with system parameters, as we will see in Sec. |Vl 



IV. DECONSTRUCTION OF INSPIRAL 
WAVEFORM GENERATION 

With the information gathered so far, let us try to un- 
derstand how the Fourier transform of the response func- 
tion is constructed, allowing for possible modifications in 
the generation of GWs. The response function is clearly 
a sum of different terms in a harmonic decomposition, 
where the term is 



2/5 

D 



(89) 



and where Q{l, 0, 0, if)) is a function of the inclination an- 
gle i and possibly all the beam-pattern functions, which 
depend on the position of the source in the sky, param- 
eterized by the angles {9, (p) and the polarization angle 
■0. The € = harmonic docs not satisfy this scaling, but 
as we have discussed in the previous sections, this mode 
contributes only subdominantly to the Fourier transform, 
and thus, it will be neglected henceforth. 

Let us now argue by symmetry and dimensional rea- 
soning that the functional form in Eq. (|89)) is unavoidable 
for the €th harmonic. If one is to have a spacetime that 
is asymptotically flat, then the metric perturbation must 
scale as D~^ to leading order. This quantity, however. 
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has dimensions of inverse length, and thus, it must be 
accompanied by a quantity with dimensions of length. 
One could use mi , m2 or a combination like the reduced 
mass fi. But since we must have that in the limit mi — >■ 
or 7712 ^ 0, then /i — > 0, one cannot normalize D with 
777,1 or 7(72, or for the same reason with Sm = mi — m2- 
Thus, the only dimension- 1 quantity available for a non- 
spinning binary is /i, which preserves exchange symme- 
try and leads to a factor of rf/^Jvi in the numerator of 
Eq. ([55)) . Also, since we will require that any modified 
gravity theory has a well-defined continuous GR limit in 
the weak field, we must disallow any other arbitrary •q'^ 
dependence. Finally, a term that scales with the har- 
monic of the orbital phase must be proportional to the 
velocity to the £*'^power. 

We should note that in principle Eq. (|5^ should be 
multiplied by an infinity series in velocity of the form 
X]ri=o ('^^," + ^t,n Inw) v" to account for higher-order FN 
corrections. In GR, one finds that for the £ = 1 mode, 
aifl = = bi Q and ai.i = = foi_i, and thus the series 
starts at IPN order with ai_2 7^ = 6i_2- In modified 
gravity theories, however, this suppression of the £ ~ 1 
mode need not be present. In fact, in theories with a 
scalar breathing mode, one finds that an £ = 1 mode is 
usually excited at leading Newtonian order. This his why 
we have not included the series dependence in Eq. ([55]) . 

Given this generic time-domain response function, we 
wish to compute the Fourier transform in the SPA. For 
this, we need both the binding energy and the rate of 
change of this energy. The former can be parameterized 
as follows 



E — Ean 



1 



A[- 

\ r 



(90) 



where we assume A is small, such that the correction 
represents a small deformation away from GR. Such a 
binding energy modifies Kepler's third law as 



1 ^ /m\P 



(91) 



or in terms of its inverse 



(92) 



Using u = roj, the Virial theorem is then modified to 



and the response function is 



(95) 



We see that this agrees with the GR and the Brans-Dicke 
results for the time-domain response function, when writ- 
ten in terms of the orbital frequency. 

The next ingredient we need is the rate of change of 
the binding energy. To model this one usually invokes 
the balance law, by which the rate of change of the bind- 
ing energy must be exactly balanced by the energy fiux 
carried away from the system (out to spatial infinity or 
into event horizons) by propagating degrees of freedom. 
Assuming that the metric perturbation accepts a multi- 
polar decomposition, the energy flux is going to be the 
sum of the square of the (^ + 1) time-derivative of the 
multipole. Let us assume a modification to the rate of 
change of the binding energy of the form 



E = Ea 



1 



r 



(96) 



where E^i-r is the GR energy flux and the second term is 
assumed small relative to the first, as we are interested 
in small deformations away from GR. Using the modified 
Kepler's law, we can rewrite this as 



E 



32 



(27rAlF)^"/^ 1 + Br]-'^""'' {2'kMF) 



29/5 



N2g/3 



-A p rr'^P/'^ {2tiMF) 



2p/3 



(97) 



We see here that there are two modifications to E: one 
coming from the modification to the Kepler law acting 
on Egji and one coming directly from the modification to 
E with the GR expression for Kepler's law. The values 
of p and q determine which of these two modifications 
dominates. We have here neglected the non-linear term 
generated by the product of terms proportional to A and 
B. 

We now have all the ingredients to compute the Fourier 
transform of h^^^ in the SPA. For this, it is convenient to 
first compute the rate of change of the orbital frequency 
via the chain rule: 



V = {mu})^/^ 



l + -Ap {mufP'"^ 
6 



(93) 



Thus, the binding energy to leading PN order becomes 
E = -]-'r]-'^''^ {2TTMFf'^ [1 



-A (5p - 6) 7?-2^^/5 {2t:MF) 



n2p/3 



(94) 
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S7y-29/5 {2TTMFf'^'^ 
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(98) 



where again we have kept terms to leading-order in the 
PN approximation and in the deformation parameters. 
The Fourier transform of the response function in the 
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SPA is then 



(99) 



where was defined in Eq. (p4| . We see then that the 
modification to the energy fiux introduces corrections to 
the amphtude of gPN order, while the modification to 
the binding energy introduces corrections of pPN order. 

This calculation shows explicitly that the leading-order 
amplitude correction induced by an £ = 1 mode will not 
be degenerate with PN amplitude corrections since the 
former enters at Newtonian order. That is, the ^ = \ 
amplitude correction in Eq. ((99|) is proportional to u-^ 
which corresponds to a — IPN order term relative to the 
GR Newtonian amplitude that scales as w"^/^. This scal- 
ing is very different from the GR amplitude corrections, 
which in Eq. (|35|) we saw scale -I-IPN or- 

der correction relative to the GR Newtonian amplitude. 
Such different frequency scalings suggest that modified 
gravity amplitude corrections induced by the £ = 1 mode 
will be weakly correlated with GR amplitude corrections. 

The Fourier phase can be similarly computed from 
Eq. (dil). We find ^'^^^ = '^'il + 5^^'\ where was 
given in Eq. ([55)1 and we have defined 



64 (4-p)(5-2p)^ ' 



15^ ^ 

64 (4-g)(5-2g) 
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(100) 



Again, the modification to the binding energy introduces 
modifications of pPN order, while the modification to the 
flux is of gPN order. 

The above expression clearly contains poles at certain 
values of q and p. These are not physical, but just a 
consequence of the assumptions made when carrying out 
the integral in Eq. ([T5|) . When p or g equal 4 or 5/2, 
the integrand that defines the Fourier phase becomes in- 
versely proportional to and thus, the result should be 
a logarithm, instead of a power law. For these special 
cases, we find 

'5*!,24,s=o = [1 - 31n^., + 5M , (101) 

'5*^to,,=4 = ^Bri-^"'4 [1 - 3 In zi, + <50o] , (102) 



'^*P=5/2,B=o = 30^77-1^,5 [1 + 3 In + 6^ , (103) 

'5*Aio,,=5/2 = y iJ^y-i^.^ [1 + 3 In ue + 5^ , (104) 

where S(t)Q is an overall constant of integration. Since 
these are 6.5PN order and 5PN order corrections relative 
to the Newtonian GR term, one usually ignores them. 



A note of caution is due at this junction. Equa- 
tion (|100[) depends only on integer powers of the reduced 
frequency ue- This is a direct consequence of our pa- 
rameterization of the energy flux and the binding energy 
in Eqs. ((96)) and (|90|) . respectively. In principle, these 
equations could be modifled by other functions that can- 
not be represented as a power series about zero velocity. 
A typical example would be logarithmic terms that, for 
example, enter at high PN order in General Relativity 
due to non-linear propagation effects (see eg. [sij and 
references therein). Moreover, certain modifled gravity 
theories could introduce screened modifications, ie. cor- 
rections that become "active" only above a certain fre- 
quency. Such effects would need to be modeled through 
a Heaviside function, which cannot be represented as a 
power-law. This, for example, is the case when dealing 
with massive Brans-Dicke gravity [33, [1^ IS] ■ 



V. PPE GENERALIZATION 

We here review the basics of standard ppE theory and 
extend this to allow for the direct presence of additional 
polarizations. We also show that for single-detector ob- 
servations, the standard ppE framework is sufficient for 
the 1 = 2 harmonic, but should be generalized to account 
for the £ = 1 harmonic component. We then generalize 
the ppE scheme to allow for multiple detectors. 



A. Standard ppE Framework 

In the standard ppE framework, one considers possible 
deformations of the two GR polarizations (h^ hx) as in- 
duced by corrections to the frequency evolution equation 
only. As we have seen, this will allow only for modifi- 
cations in the £ = 2 harmonic of the Fourier transform. 
Such corrections can arise due to deformations to the 
binding energy of the binary or the energy flux carried by 
all degrees of freedom away from the binary. In Sec. IIV| 
we computed such corrections to leading PN order and to 
leading-order in the deformation parameters. As Yunes 
and Pretorius found [l^ . the corrections to the Fourier 
transform of the response function can be generically cast 
as 

Ve,o(/) - A U2-'/'e-^*Gl(i + au«)e'2/3"^ , (105) 

where we have inserted a factor of 2 in the phase cor- 
rection for future convenience. If the modiflcations to 
the binding energy enter at the same PN order as the 
modiflcations to the energy flux, ie. q = k = p, then 



{6 + 4k- 5fc2) - is 



-2k/5 
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64 (4-fc)(5-2A:) 
a = 2k, b = 2k-5. 



ri (106) 

[A (5fc2 - 2fc - 6) + 35] 77-2^/5 ^ 

(107) 
(108) 
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If p < q, then to leading-order one recovers again 
Eq. PI13|) with 



6 



(109) 



P = -^7T^-^^^A (5/ - 2p - 6) ,r'^/' : 



64 (4-p)(5-2p)' 

a = 2p , b = 2p — 5 . 
while if g < p, then you have Eq. (|105p with 
1 



(110) 
(111) 



2 ' 



/3 



15 



(112) 

B?7-2«/5 . (113) 



64 (4- g)(5 - 2q)' 
a = 2(7, 6 = 2g-5. (114) 

We see then clearly that all possibilities can be mapped 
to the ppE waveform family of Eq. pOSp . 

The ppE waveform family of Eq. (|105p depends on the 

standard 5 system parameters Agr = {A,A4,ri,tc,^c ) 
and on 4 ppE theory parameters AppE,o = (a, O; b). The 
ppE parameters (a, 6) are pure numbers that take discrete 
values and characterize the type of physical modification 
to the GR prediction. The parameters (a, /3) depend on 
the coupling constants of the modified gravity theory, as 
well possibly on the masses and spins of the binary, and 
they are the quantities we wish to constrain with observa- 
tions. When (a, f3) = 0, one recovers GR, while for other 
values of ppE parameters, one recovers the predictions of 
other theories for the £ = 2 harmonic. For example, the 
Brans-Dicke prediction for the £ = 2 harmonic is recov- 
ered when 

224 5 

«Bo = ^/3bo , /3bo = -^V'^'S^^ , (115) 
3 7168 



^BD = -7. 



(116) 



Rosen's and Lightman-Lec theory cannot be easily 
mapped because they do not possess a GR limit, and 
thus, they are not technically small deformations away 
from Einstein's theory (see Sec. 

From a physical standpoint, Eq. (|105p makes perfect 
sense: the waveform depends on an amplitude and a 
phase degree of freedom, so any deformation of it can be 
mapped to a change in these two quantities. One could 
in principle Taylor expand these functional degrees of 
freedom, but lacking a particular theoretical framework, 
the controlling factor of such an expansion is unknown. 
Thus, since different theories predict different controlling 
factors, it is best to use the parameterization in Eq. (|105p . 
The pure numbers (a, 6), however, should not be thought 
of as completely free, as after all they are generated by 
velocity corrections to GR's physical principles. If one as- 
sumes that modified gravity theories always lead to cor- 
rections of the form v", where n E 1j, then (a, 6) G Z. 
To date, we have not encountered any modified gravity 



theory that predicts corrections in non-integer powers of 
the orbital velocity. 

Can modifications to the radiation-reaction force or 
the Hamiltonian in modified gravity theories depend on 
non-integer powers of the velocity? We believe the an- 
swer to this question is no. If such terms were present, 
the corrections would be non-analytic in the weak-field, 
ie. they would lack a well-defined Taylor expansion about 
zero velocity. Such lack of analyticity would lead to a 
breakdown of the initial value formulation. For exam- 
ple, differential equations of the form v = with p a 
non-integer and initial condition z;(0) = do not possess 
a unique solution. Notice that this is not a problem for 
terms proportional to negative powers of velocity or to 
natural logarithms, as the latter are analytic everywhere 
in their domain (0, oo), ie. f = is not in the domain of 
these functions. 

However, it is known that non integer powers of the 
frequency in the phase of the GWs arise when the orbit 
is no longer circular. The reason the mathematical argu- 
ment presented above is no longer valid is the following. 
For eccentric orbits, the binding energy is a bivariate se- 
ries both in the azimuthal velocity U0 and in the radial 
velocity Ur- The former is initially non-zero, but very 
close to zero, increasing as the inspiral proceeds. There- 
fore, one can expand the binding energy around ug = 
and a non-integer power in the evolution law would break 
uniqueness. On the other hand, during the inspiral Ur de- 
cays. An inspiral can begin with large eccentricity, but 
due to radiation reaction, it would lose eccentricity and 
become circular as the inspiral proceeds. Therefore, one 
can no longer expand the binding energy around Ur = 0, 
which implies the initial condition Ur{0) = is no longer 
valid. A differential equation iir — with p non-integer 
and an initial condition Ur(O) 7^ has no uniqueness 
problem. Therefore, if an extra field (j) were present and 
it decayed during the inspiral without an initial condition 
0(0) = 0, one cannot use the argument above to rule out 
non-integer powers in the evolution of the GW phase. 

The parameterization of Eq. (|105p does not use the 
fact that a modification to GW generation will usually 
introduce changes to the amplitude and phase that will 
be related to each other. Indeed, from the analysis of 
Sec. IIV) we see that a parameterization of the i = 2 har- 
monic of the Fourier response that explicitly recognizes 
these relations would be 
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Notice that a and a have been completely eliminated. 
The parameter c S M is not a free parameter, but it is 
fully determined by b. When a conservative correction 
dominates, then 



Cnnns 



16 6(3-6)(42& + 61 + 5&2 



15 562 + 46& + 81 
when a dissipativc correction dominates, then 

16, 
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and if both conservative and dissipative corrections enter 
at the same PN order, then 



Cboth — 



32 5(3-6)(446 + 71 + 552) 
15 562 + 465 + 81 



(120) 



Therefore, the ppE theory parameters of /ippE,i are 

^ppE.i = (fe, two less than AppE,o- This might seem 
strange since originahy we started with 4 arbitrary con- 
stants {A, B,p, q), two to parameterize corrections to the 
binding energy and two for modifications to the flux. 
Without loss of generality, however, either one of these 
two will be dominant or they will enter at the same PN 
order, enabling us to eliminate 2 parameters. Notice that 
when b = —7 and a dissipative correction dominates, then 
Cdiss = 224/3 and one recovers the Brans-Dicke result of 
Eq. dm]). 

This reduced parameterization, however, neglects pos- 
sible modifications in the propagation of GWs, due for 
example to a Lorentz- violating graviton dispersion rela- 
tion. The analysis of [4l|, \^ showed that modifications 
in GW propagation only introduce corrections to the 
Fourier phase, and thus, Eq. (|117p should be modified 
to 



(121) 

where k E Z while k G M depends on the coupling 
constants of the theory. Obviously, if one neglects the 
2/3^2 term in the phase, one recovers exactly /ippE.o 

af- 
ter an appropriate relabeling of the ppE parameters [see 
Eq. pUS)) ]. This is only vahd if fc < &, as then the 2^w^ 
term would be of higher PN order than the term. 
This template family has parameters AppB,2 = (6, k, fc), 
the same number as AppE,Oj although there is a clear phase 
degeneracy here when k = b. The ppE template family 
/ippE,2 allows one to disentangle generation from propa- 
gation effects. 

The paramcterizations discussed above, however, only 
work for a single-detector and for the £ ~ 2 harmonic. 
In the next subsections, we generalize these proposals 
to allow for power in the £ = 1 harmonic and for the 
possibility of multiple-detector detections. 



B. Generalized ppE Scheme: 
Single Detector, all harmonics 

Given a single interferometric detector, one cannot 
break the degeneracies in the waveform to separately ex- 
tract all physical parameters [l^ . For example, one can- 
not separately measure all the angles that characterize 
the position of the source in the sky. Thus, even within 
GR, the prototypical template family is a reduction of 
Eq. (1221) to 



-7/6 g-i'I'GB 



(122) 



where we have absorbed a constant factor into the phase 
of coalescence. That is, one can measure the parame- 
ters contained in ^'gr, ie. {tc, $c,?7,-M), plus the overall 
amplitude Agr, which essentially determines the signal- 
to-noise ratio. However, one cannot separately measure 
all the quantities encoded in 



Ac 



/57r^ 




V96, 


' D . 



Fl{l 



cos^ l) 



AF. 



? 9 

cos L 



1/2 



(123) 

Similarly, when deciding how to modify the GR tem- 
plate family to allow for modified gravity effects, one 
must keep possible degeneracies in mind. Below, we cal- 
culate the reductions of the SPA Fourier transform of the 
response functions of Sec. IIIII for each modified gravity 
theory studied and then discuss a few ppE proposals. 



1. Brans-Dicke Theory 



Equation (j53p can be rewritten as 



^SD (/) — A^T>U. 



-7/2 
2 



. 224^ 

1 + — PbdW2 ■ 



-9/2 - 

-7bdWi e 



(124) 



where ^-^^i^ = + S^^^^l is given by Eq. with the 
substitution i>c -> <i>c in ^'gr- The Brans-Dicke correc- 
tion can be written as J^'bd = ^ (^/2)/3bdm7''' where /3bd 
was already defined in Eq. (jllSp and we have redefined 
the chirp mass M = M[l - (^/15)(rVl2 - k^^)]. In 
rewriting the waveform in Eq. (j53p in terms of an over- 
all amplitude and an overall phase, we also redefined the 
(constant) phase of coalescence as 



3^,2 



arctan 



2 cos iFy 



(1 + cos2 l)F+ 

^rcost(i -cos^ QFxJ; 

' (1 + cos2 i)2F2 + 4 cos2 lF^ 



(125) 



which leads to two distinct constants for each harmonic. 
Therefore, the parameters of such a template family are 
the usual system parameters Agr, plus the Brans-Dicke 

parameters Abd = (/3bd, 7bd, *i^^)- 

For a single-detector, one still detects only a combined 
overall amplitude ^bd for the £ = 2 harmonic, which in 
terms of the fundamental parameters of the theory and 
the bcani-pattcrn functions is given by 



-^BD 



Fld + cos^ i)^ 



AFI cos^ l 



F+F^ {cos* i - 1) e r] 



1/2 



(126) 



and we recall that F is defined below Eq. (j4T|) . We have 
here absorbed a prefactor into the luminosity distance, 
namely D = D [l + ^ (F2/24 + A:bd/2)] , where ^bd and 
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S are defined in and around Eq. (|4T|) . Moreover, for a 
single-detector one cannot measure ^ separately, but the 
combined phase and amplitude quantities. The latter is 
related to ^ via 



7bd 



48 



D 



i^/^S^F^ Sim, 



(127) 



Wc separate 7bd from /3bd here because they appear mul- 
tiplied by different powers of Uf and the beam-pattern 
functions in the Fourier transform of the response func- 
tion. That is, different (single) detectors will measure 
the same /3bd, but different "/bv- 



two sets and theory parameters is 



Q [Fl + {F^ + Fl) sin^ 



+ Fl cos" L + 2F,Fl sin^ l+{F^ + F^) F^^ sin 2l] ^'^ , 

(131) 



(132) 



As in the Brans-Dicke case, all detectors will see the same 
/3r for a given event, but different 7^. 



2. Rosen's Theory 



Lightman-Lee Theory 



Equation ()69p can be rewritten as 



where ^i^^ is given by Eq. (|72p with the substitution 
^•c '^c ^- The last term in Eq. ^ can also be written 
as where /3r is given below. As in the Brans- 

Dicke case, the above expression is obtained by rewrit- 
ing Eq. in terms of an overall amplitude and phase, 
where we have redefined the (constant) phase of coales- 
cence via 



■7rU. 



(128) 



Equation 



can be rewritten as 



-7/2^-1* ■ 



+ 7llUi 



-9/2 



(133) 



where ^'[l'' in this equation is given by Eq. ((88)) with the 
substitution $c ^ The last term in Eq. dHU) can 

also be written as where /3ll is given below. In 

rewriting Eq. (j85p in terms of an overall amplitude and 
phase, we had to redefine the constant phase of coalesce 
to 



$(^) = <l)^ — (5^^ arctan ( — 



Fb sin L + F, sin l + F„. cos l 



F„ 



• 2 arctan 



i^+(l + cos2 i) 
2Fx cos L — 2F,„ sin t 



(F, +F0sin2i-HF3,sin2i 



2Fx cos L — 2K„ sin l 



(129) 



Therefore, the parameters of this template family are the 
usual system parameters Aqr, plus the Rosen parameters 
(/?,,7b,$W). 

As before, a single detector can only measure the com- 
bined amplitude ^r, which is related to other system 
parameters via 



4:D 



3/4 r 



Fi(l + cos^ if + AFl cos^ i 



+ Fb sin'^ L + AFl sin'' l + F^ sin^ 2t + AF^ sin^ t 

+ 8 cos i sin^ + 4 sin'' tFtFL 4- 4 cos t sin^ tFbF„ 

-f 2(l-cos4 l) (F+F, + 2F+F,) 

+ (1/2) (6sin2t-|-sin4t)F+F,, +4sin2tFxF„//^ , 

(130) 

One can also measure the combination of theory con- 
stants in /3r and 7r, where the relationship between these 



= $c - 5i,2 arctan 



F4- (1 



- cos^ l) 



2F,„ sin L + 2Fx cos l 



(2Fl - 3FJ sin^ l + F,, sin 2t 



Si 1 arctan 



2F,„ sin t -I- 2Fx cos t 

X5/4)Fb + Fl] sin t - cos lF,, 



F.. 



(134) 



Therefore, the parameters of this waveform arc the usual 
ones Agr, plus the Lightman-Lee ones (/3llj 7llj 'i'c^^)- 

As before, a single detector can only measure the over- 
all amplitude v4ll, which is related to other parameters 



'''' ^ V 168~D~ 



Fl (1 -Fcos^ t)%4cos2 lFI 



- (9F2 + AFl) sin^ L + sin^ 2lFI + 4 sin^ lF^ 

- 8 cos t sin^ iF^F,, - 12 sin'' lF^F^ - 12 sin^ t cos tF^F,, 

- (1/2) (5 sin L + sin 3i) F+ (3 sin iFt, - 2 sin tF^ 

-2cosiF,<,) + 4sin2iF3„Fx]^^^ , (135) 



Moreover, one can also only measure (7ll, /^ll), and these 
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are related to other fundamental theory parameters via 

'25 



7l 



189 ' D 



+ F^ ] sm^ , 



cos^ tF^ + Fl + - sin^ lF^F^ 



-.1/2 



/3. 



625 
16464 



4 



F^F^^ sin 2t 



(136) 
(137) 



Multiple detectors will measure different 7ll, but the 
same /3ll given the same source. 



4- ppE Scheme 

The above discussion allows us to conclude that, al- 
though the already existing ppE scheme is perfectly ade- 
quate to capture the 1 — 2 harmonic in the Fourier trans- 
form, it docs not capture the ^ = 1 harmonic. Power in 
this harmonic is a natural consequence of theories that 
predict extra polarization states. The largest general- 
ization of the ppE scheme that covers all alternatives is 
simply 



Ce.o(/) 



+ 7 iti e 



(138) 



where now the theory parameters arc AppE,o — 

(a, a, /3, 6, 7, c, (5, d, $c^^). We say that this is the largest 
generalization because, in addition to the standard 5 sys- 
tem parameters Aqr, we have added 9 more for a total of 
14 parameters, without imposing any relations between 
them. 

One might worry that the parameterization of the 
1=1 term (the second line in Eq. ()138|) will be degen- 
erate with PN amplitude corrections. This will indeed 
be the case if c e (—5/2, —2, —3/2, . . .). Modified gravity 
theories that predict the excitation of a scalar breathing 
mode, however, are likely to introduce £ = 1 corrections 
to the waveform that enter at leading, Newtonian order. 
For these, c = —9/2, which is a — IPN order correction 
in the amplitude relative to the leading-order GR term. 
As such, these modifications should be weakly correlated 
with PN amplitude corrections. 

Since the leading-order term of the I —1 harmonic en- 
ters at Newtonian, leading order, the amplitude parame- 
ter 7 has to be very small, for the theory to be consistent 
with binary pulsar observations. Given that detectors 
are more sensitive to the GW phase than its amplitude, 
it would appear that the effect of additional polariza- 
tion states would first show up indirectly in the modified 
phase evolution of h+ and , and not through the direct 
detection of the additional polarization states. However, 
E ^ 0{t^), where e is the small parameter mea- 

suring the deviation of the theory in question from GR. 



Then, it is clear from the SPA formalism (Eq. (iMl) ) that 
the amplitude of the Fourier transform of the i = \ har- 
monic will have 2 terms: one coming from the quadrupole 
part of F and one coming from the dipole part of F. The 
former is of order e, while the latter of order and can, 
thus, be neglected. From Eq. ([15]) we see that (3 is propor- 
tional to and we still have the possibility of measuring 
7 before /?. Only in very specific theories, for example 
in Brans-Dicke theory where the kinetic scalar term is 
multiplied by the reciprocal of the Brans-Dicke coupling 
parameter, is E modified from its GR expectation at lin- 
ear order in the coupling. A more detailed study of the 
possibility of measuring 7 is left for future work. 

In the above parameterization, we did not use the fact 
that some of these parameters are not independent, if 
they are all to arise from a modification to the generation 
of GWs. As we showed in Sec. IIVI given a correction to 
the energy flux or the binding energy that scales as 
relative to the leading-order GR expectation, the correc- 
tion to the phase in the ^ = 2 harmonic is proportional 
while the correction to the amplitude goes as 



to u 

„,2k 



2fc-5 



One can then postulate a reduced ppE model of the 



form 



Ce,i(/)=-4«2 



7/2 



(2) 



[l-f c I3u. 



b+5j gi20«° 



(139) 



where is given in Eq. with $c -> We have 
also explicitly used the fact that the f = 2 harmonic 
leads to an amplitude proportional to /~^/®, while an 
£ = 1 harmonic leads to f~^^'^ when only the leading- 
order quadrupole radiation is taken into account. No- 
tice that c is again given by Eqs. (jllSp - (jl20p . ie. it 
is fully determined by the conservative or dissipative 
modifications. The ppE theory parameters of h^^^ 1 are 

AppE.i = /3, 7, *&c^''), in addition to the usual system 
parameters Aqr, for a total of 9 parameters. 

The result in Eq. (|139p is very similar to that of 
Arun 24 1, so let us compare it directly. In [23], Arun 
considered the effect of dipole radiation in the Fourier 
waveform, which corresponds to our waveform with b = 
—7. Moreover, he neglected any conservative corrections 
(modifications to the binding energy), and thus, c = Cdiss 
as given in Eq. (|119p . In this limit, h^°^ ^ reduces to his 
Eq. (9). We note that Arun's overall amplitude is com- 
plex, and thus, there is an additional parameter that was 
not made explicit in his equations. We have here pulled 
this factor out by including ^c^' ^ f&c^-' explicitly. 

Whether one chooses Eq. (|138p or (|139p depends on 
the question that one is trying to answer. The param- 
eterization in Eq. (|139p contains fewer parameters, but 
stronger assumptions about the type of corrections one 
is searching for. In deriving this equation, we have ne- 
glected all corrections that cannot be represented by a 
power-law. For example, we have ignored possible loga- 
rithmic modified gravity corrections in the Hamiltonian 
or radiation-reaction force, as well as screened modifi- 
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cations that depend on the Heaviside function. More- 
over, we have also neglected the fact that GR modifica- 
tions consist of an infinite power series in velocity, and 
thus, corrections to the phase will not depend on a sin- 
gle term, but an entire PN series. The parameterization 
of Eq. (|139p attempts to extract only the leading-order, 
controlling factor in the expansion. If the true theory of 
nature contains these type of deviations, Eq. (|138p might 
be more powerful at detecting a deviation. This is the 
topic of a follow-up study that is in preparation [4^ . 

Finally, we can also allow for propagation effects by 
multiplying the phases in Eq. (|139p by the appropriate 
ppE propagation factors: 



2if)=Au, 



-7/2 



(2) 



"l + c I3u° 



I -9/2 - 



-,(1) 



(140) 



where (k, k) are new ppE propagation parameters. Us- 
ing such a waveform, however, might not be ideal for 
data analysis purposes, as usually either u'' or u'^ will 
dominate, and it might be preferable to ignore the sub- 
dominant contribution. 



C. Generalized ppE Scheme: 



Multiple Detectors, all harmonics 

The theory amplitude parameters, such as /? and 7, in- 
troduced in the previous subsection, clearly depend both 
on the coupling constants of the theory, as well as on 
the symmetric mass ratio, the sky location of the source 
and the inclination angle. Therefore, these parameters 
are system and detector dependent and not suitable for 
comparing data acquired from a detector network. In this 
subsection, we will construct a generalization to the ppE 
parameterization for the waveform amplitude to resolve 
this problem. 

Let us then begin by presenting the Fourier transform 
of the response function in the SPA for all the theories 
considered in Sec. lIIIl but in full detail this time, without 
rewriting the response in terms of an overall amplitude 
and phase correction. 

1. GR : 



hZif) = -[^^+(1 + cos^ + 2i^^x cost] 



96 



1/2 



M 



D 



(141) 



2. Brans- Dicke Theory: 

h^M) = (-[^^+(1 + cos^ i) + 2iF^ cost] 



p2 _|_ ^BD 



e s' 



[F+{1 + cos^ l) + 2iF^ cosi] y^}n 



-7/2 



Fb sin t - ^ ^2 



7/2 



96 



D 



e BD 



(142) 



3. Rosen's Theory: 

KM) = {- [^+(1 + COS^ + 2»Fx COS i] 

+ [-Ft, sin^ L - 2i^L sin^ l - 2F,J sin t - F,, sin 2i\ } 



336 



1/2 



-3/4 



D 



^ (-F, sin. - F, sin. - F.„z - F.. cos.) gk'' ^^^'/'^u-"\-^'<' , 



(143) 



4. Lightman-Lee Theory: 

Kin if)^ + cos^ .) + 2iFx cos l] + [3Fb sin^ . - 2F^ sin^ l - 2F,J sin . - sin 2l] } 

"5 



577 M'^ -7/2 _,>]/(2) 



336 



D 



-Fb sin i — Fi^ sin l + Fg„ + F^^ cos . 



(144) 



where the phase terms are given in Sec. 

Using the above examples for guidance, one can generalize the ppE framework to apply to multiple detectors via 
the large ppE class 
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+ [7+F+ + 7x i^x + Tbi^b + 7l^L + %„i^s„ + 7?i/^«^e-'*o^e'*< . (145) 



Comparing this to Eq. (|138|) . it is obvious that the pa- 
rameters (a, /3, b, c, S, d) are the same as introduced for 
a single-detector. In order to account for muhiple de- 
tectors, however, we have introduced a number of pa- 
rameters in place of (a, 7, one for each polar- 
ization mode: a — > (q!-|_, ax , Q^b, Q^l, 0^50, ttan) and 7 — > 
(7+, 7x , 7b, 7l, Tao, Tan)- Wc have here absorbed any pos- 
sible inclination dependence into these parameters, as 
multiple detectors will essentially see the same inclina- 
tion angles. Also, we have here included both a+ and 
a X for the sake of generality, but in all examples encoun- 
tered usually these parameters are not independent. 

The parameters introduced are either purely real or 
purely imaginary, but not complex. Wc can then rewrite 



them as a, = (5^6*'^/^" and 7^ = 7^6*'^/^^ where {ai,^i) 
are real parameters, while (di, 7^) are binary parameters, 
ie. they are either equal to unity or they vanish. Although 
this introduces a new parameter d or 7 for each a or 7 
parameter, the former two are not full parameters, in that 
they only increase the volume of the parameter space but 
not its dimensionality. 

As we discussed in Sec. lVBI and llVi however, not all of 
these parameters are independent. Using the arguments 
of Sec. IIVI we can restrict the large class of templates 
of Eq. (|145[) to the leading-order corrections induced by 
modified GW generation corrections without preferred 
frames: 



K"^M) = h'Zil + c f3v!'+^)e^'^< + Ki^bSin^.-f a,i^,sin2t + a.„F.„sini + a3.F.„sin2t]:^7.2-'/2e-'*ole2*''"^ 



+ [^,F, sin L + 7,F, sin l + x„F.„ + 7..F.„ cos i] r^^/^^u^^'^^e-'^oke'P^ . 



(146) 



As before, there is no need to introduce $c 
here. The 10 theory parameters are then 

(&,/3, ab,aL,aao,a=,„,7b,7L,73c,7s„)- Of these, 6 e Z, 
while all others depend on the fundamental coupling 
constant of the theory. 

We have proposed the specific parameterization of the 
amplitude corrections based on the following considera- 
tions. For theories with no preferred directions and for 
circular orbits, the leading-order, trace-reversed metric 
perturbation can have terms proportional to 

1. For^ = 2 

- {N ■ vf , {N ■ xf (147) 

- {N ■ v)v' , {N ■ v)x^ , {N ■ x)v' , (TV • x)x' (148) 

- {N ■ vfS'^ , {N ■ xfS'^ , v'v^ , x'x^ (149) 

2. For ^ = 1 



If we apply the projectors in Eq. @ to the above equa- 
tions to extract the polarization amplitudes and then use 
Eq. ^ to calculate each mode, we find that the incli- 
nation dependence given in Ea. (|146p is unique for the 
leading-order terms. This is confirmed in the theories we 
discussed above, all of which lead to the same structure: 
the £ = I harmonic has a term proportional to Fa and 
a similar term appears in the £ = 2 harmonic but multi- 
plied by another factor of sin t and a different amplitude, 
where the latter depends on the particular theory con- 
sidered. The inclination dependence will be different for 
terms that enter at higher order in the orbital velocity, 
just as happens in GR. We have here not allowed for 
ppE theory parameters multiplying F+ and Fx in the 
£ — I harmonic because such terms could only arise if 
hij ^ Xi Aj + Vi Bj , where {Ai , Bi ) are independent of 
the orbital phase, which would be indicative of theories 
with preferred directions. Such theories exist [TT| - [l3| . 
but we leave a detailed analysis of these effects to future 
work. 



h"" - {N ■ v) ,{N-x) (150) 

h^i X-' (151) Finally, we can also allow for propagation modifica- 

_ . . ^ ^ . . .. tions, in addition to the generation ones, by introducing 

h ^ ^ (TV • v)S ^ , (TV • x)S (152) ,^ further phase corrections of the form 
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+ [ofbi^b sin'^ L + Q^L^L sin^ l + ia^^F^„ sin l + a^^F^^ sin 2t 



+ [7bFb sin t + 7lFl sin t + i7,„F,„ + %,F,, cos i] r^i/^^u^ ^/^e-'^^ol 



(153) 



where (k, k) are propagation ppE parameters. As in 
Eq. (|140p . using this waveform might not be ideal, as 
usuaUy either or u*"' will dominate, and one might wish 
to neglect the subdominant contribution. 

Before proceeding, let us make a remark about the GR 
limit of the parameterizations discussed here. We saw 
earlier that both Rosen's and Lightman-Lee theory do 
not possess a continuous GR limit. That is, if one takes 
the coupling parameters of this theory to zero ^ — ^ 0, 
one does not recover Einstein's theory. While our pa- 
rameterization is designed to have a smooth GR limit, 
it does manage to cover these exotic and observationally 
disfavoured alternatives. More precisely, the parameteri- 
zation discussed above poses no restriction on the magni- 
tude of the amplitude /3, and thus, if /? were large enough 
and b ~ —5, the phase terms proportional to /3 could 
combine with the GR terms to change the leading-order 
behavior. Such behavior can be discarded if one imposes 
priors on /3 and k, for example, by requirin g th em to 
satisfy their current binary pulsar constraints j44| . 



VI. NULL STREAMS 

A general method used in GW data analysis to sep- 
arate signals from noise is that of null streams. This 
approach was first introduced by Giirsel and Tinto [45| 
and was later extended by Chatterji et al. [4^. The idea 
behind this is to combine the data from a network of de- 
tectors to find linear combinations that contain no GW 
signal, only noise, in the hopes of separating false-alarms 
from real GW burst events. Their analysis only included 
the /i+ and hx polarizations, since they were considering 
only GR events. 

We can extend this method to test Einstein's the- 
ory through the proper combination of output from sev- 
eral detectors. In essence, one can extend the analysis 
of to allow for all possible polarizations and then 

construct the appropriate null streams ( "null" within 
GR) through appropriate projections. This method is 
promising to test GR, since the detection of a GW in 
a GR null stream would automatically signal a devia- 
tion from GR (assuming the source location is known 
- the more general analysis for sources with unknown 
sky location will be developed elsewhere). One can thus 
search for statistically significant deviations from noise 
in GR null streams, both with a template (as given in 
Sec.|V| and without one. Given a detection of a signal in 



a GR null stream, one could then reconstruct the signal 
through the templates in Sec. |Vl 

Let us assume that there exist D > 6 detectors with 
uncorrclated noise and that, for a given source, we know 
its position in the sky, as might be the case if we have 
an electromagnetic counterpart. Given this, one knows 
exactly how to time-shift the signal from detector to de- 
tector. For a detector a, the noise- weighted signal from a 
source at location Cls in the sky in the frequency domain 
is 

da = F+h+ + Ffhx +FfK^ + Fak. + FfK + FfK + fia ■ 

(154) 

Here, and in the remaining section, d^, F^ and fia are 
the noise-weighted signal, the antenna patterns and the 
noise of the ath detector respectively, each defined as 
their standard value divided by y/SJJ)/2 where Sa{f) 
is the power spectral density. 

Given D detectors, we can then rewrite Eq. (|154p as 
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or using tensor notation. 



(155) 



(156) 



where the a runs over the number of detectors and the j 
over the polarizations. The quantity F"" j is acting anal- 
ogous to a metric in the signal manifold, and thus, the 
first term on the right-hand side can be interpreted as 
the projection of the wave vector W along the directions 
F+a,F''a,F'\,F'-a,F\ and F^. Therefore, one can 
create data sets that have no component of a certain po- 
larization by projecting them to a direction orthogonal 
to the direction defined by the beam pattern functions of 
this polarization mode. This is illustrated in Fig [T] for 3 
detectors. 

For D detectors, the signal manifold is D dimensional 
with D basis vectors. 5 of which can be chosen along 
the F+a, F'^a, F'^a, -F'"a and F^'a directions. The reason 
why we cannot choose 6 linearly independent basis vec- 
tors along the F^ directions is that F^ = —F^, as can be 
clearly seen from Eq. [2] (see also [13, Ell). This means 
that although we have a 6 x D matrix, this only has 5 
linearly independent columns. The remaining D — 5 vec- 
tors would give us complete null streams, ie. streams with 
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FIG. 1: _ Schematic diagram of the projection of the data 
stream dc orthogonal to the GR subspace {F^ , F^^ ) for 3 de- 
tectors to build the GR null stream dcR.nuii- 



truly no signals, no matter how many polarizations arc 
present. Currently, there are only 3 detectors active or 
soon to be active and collecting data, and thus, we can 
only eliminate two polarization modes from the full sig- 
nal. Let us then define the GR null stream as the stream 
that has no tensor modes {h^ and hx), namely 
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(this is analogous to the null stream defined in |46[). 
Clearly, this will contain a mixture of the two scalar 
modes and two vector modes. 

Given data from 3 detectors, one can then study 
whether (iGR,nuii contains any statistically significant de- 
viations from noise without ever invoking a template. 
However, since we have constructed a parameterization 
for a response function with all polarizations, one could 
filter 

•^GR.nuU with the GR null scalar 



rir 
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fitting for the system parameters, together with the ppE 
theory parameters. The quantity he here should be in- 
terpreted as, for example, /i™ i given in Eq. ()145|) for 
the cth detector. 

What if we had more than 3 detectors? Each addi- 
tional detector would allow us to eliminate one of the two 
vector modes or both the scalar modes, beyond the two 
tensor ones. For example, given four detectors, (icR null 
becomes a two vector, as there are two GR null streams. 



which in addition also contain no power in one of the 
vector modes or the scalar modes, depending on which 
direction we chose to project out. Therefore, given 5 de- 
tectors one can construct 3 enhanced GR null streams, 
each with power in a signal direction. With more than 5 
detectors, one can then construct complete null streams 
(in any theory of gravity). This topic and its implemen- 
tation will be studied more carefully elsewhere [43j]. 



VII. CONCLUSIONS 

We have studied how theories that predict the exis- 
tence of all GW polarization modes affect the GW re- 
sponse function. First, we considered three modified the- 
ories (Brans-Dicke theory, Rosen's theory and Lightman- 
Lee theory) and, for each of them, we extracted the po- 
larization modes and calculated the response function 
in the time and frequency domains, for a quasi-circular, 
non-spinning compact object inspiral source. We found 
that, although generically the response function contains 
terms proportional to the £ = (0, 1,2) harmonic of the 
orbital phase, only the i = (1,2) harmonics contribute 
significantly to the Fourier transform in the SPA. This 
is because the £ = harmonic lacks a stationary point 
in the generalized Fourier integral, and thus, it becomes 
subdominant relative to the £ = (1,2) harmonics by the 
Riemann-Lebesgue lemma [2l| . 

We then considered the generic structure of the £- 
harmonic response function in the time domain using 
symmetry arguments and dimensional analysis. We 
Fourier-transformed this generic structure in the SPA 
and found that the leading-order Fourier amplitude and 
phase scale as £~^l'^rf'^~^^l^%Li^ '^'^)/^ g^j^^j respec- 
tively, where ui — (27rA^//£)^/'^. 77 is the symmet- 
ric mass, / is the GW frequency and Ai is the chirp 
mass. We then allowed for power-law modifications to the 
radiation-reaction force (the energy fiux) and the Hamil- 
tonian (the binding energy) . Given a modification of rel- 
ative nth PN order to either of these quantities, we found 
that the Fourier amplitude and phase acquire corrections 
that both scale as u^" relative to the GR leading-order 
term. 

With this in hand, we considered the ppE framework 
and whether it was capable of capturing this higher- 
harmonic modifications to the response function. We 
found that, although the original ppE framework is per- 
fectly capable of modeling the £ = 2 harmonic, it needs to 
be modified to also allow for the £ = 1 harmonic. We pro- 
posed two possible extensions to the single-detector ppE 
framework: a generic one with a total of 9 ppE theory 
parameters; and a reduced one with only 4 ppE theory 
parameters that uses some of the scaling relations found 
for a generic terms promotional to the ^-harmonic of the 
orbital phase in the response function. 

The original ppE scheme was not applicable to multiple 
detector configurations, as it was based on a generic de- 
formation of a GR single- detector response. When trying 
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to constrain the existence of non-GR polarizations, how- 
ever, one needs to redo this analysis to allow for multiple 
detectors. We carried out such an extension proposing 
another two ppE, multiple-detector template families: a 
generic one with a total of 18 ppE theory parameters; 
and a reduced one with a total of 10 ppE parameters. 
The large increase in ppE parameters is necessary if one 
wishes to test for the existence of all 4 additional non- 
GR polarizations in both the £ = 1 and £ = 2 harmonic 
of the Fourier transformed response function. If one re- 
stricts attention to only one additional non-GR polariza- 
tion, then one requires only 9 ppE parameters for the 
generic proposal and 5 ppE parameters for the restricted 
proposal. 

Some of the parameters we introduced here are already 
constrained by binary pulsar observations. For example, 
if there exists a scalar breathing mode at leading, New- 
tonian order, then it will modify the amplitude through 
a — IPN order term, which then will change the rate of 
change of the binding energy, and thus, of the orbital pe- 
riod. Such negative-PN order terms must be very small 
in magnitude to be compatible with binary pulsar obser- 
vations 441. Such restrictions could be accounted for in 



the form of a prior through a Baycsian analysis. 

We conclude with a discussion of a possible method to 
test for the existence of additional non-GR polarization 
through the combination of the output of multiple detec- 
tors. That is, we extend the concept of null streams to 
non-GR theories, defining a GR null stream as one that 
would be consistent with noise in GR, but that would 
contain non-GR polarizations in a modified gravity the- 
ory. Such a null stream is simply constructed by build- 
ing basis vectors orthogonal to the 2 GR polarization 
basis vectors in the signal manifold, and then project- 
ing out the response with these orthogonal directions. 
We find that at least 3 detectors are necessary to con- 
struct one such null stream. One requires 5 detectors to 
isolate the 4 additional non-GR polarizations into single 
GR null streams and 6 detectors arc required to create a 
null streams in all theories of gravity. 

Although the analysis presented here is a definite step 
forward toward the construction of model-independent 
tests of GR with GWs from binary inspirals, it leaves 
open several directions for future research. The four ppE 
proposals discussed in this paper were in part inspired 
by our analysis of three specific models. In principle, 
one would like to repeat the analysis carried out here 
for vector-tensor theories, such as @, stratified theo- 
ries H, [l3- 1£ I and tensor- vector-scalar theories, such as 



TeVeS [9|, [lOf or Einstein- Aether theory [llHl3|. The 
latter will be particularly helpful in improving the ppE 
scheme, as it predicts the existence of preferred direc- 



tions, which should introduce £ = 1 harmonic modifica- 
tions to the time-domain response function. Another step 
towards the improvement of the above proposals would 
be to determine a generic parametrization for the quan- 
tities that depend on the parameters of the system, like 
the mass, the spins etc. 

One of the main motivations for the ppE scheme is to 
develop a framework to test GR in the strong-field and 
here we have concentrated on the inspiral phase only, 
which might seem like a contradiction. As explained 
in the Introduction, however, by "strong-field" we here 
mean the region of spacetime where the gravitational 
field is dynamical and non-linear. The PN expansion 
allows one to model such non-linearities and strong dy- 
namics perturbativcly, and it is precisely this that the 
ppE framework modifies. Ideally, one would also want 
to consider extensions of the ppE framework during the 
plunge and merger phase. In these phases, one expects 
the most amount of radiation emitted through non-GR 
polarizations, and perhaps, the largest deviations from 
our GR expectations. The plunge and merger, however, 
are extremely difficult to model even within GR. Lacking 
a numerical modeling of these phases in modified gravity 
theories, one can only improve on inspiral ppE models 
at the current time. In the future, however, it would 
be most interesting to see how additional polarizations 
modify GR waveforms during merger. 

Another obvious direction for future research is the im- 
plementation of such proposals in a realistic data analysis 
pipeline, like that discussed in [l^ or [i^. In particu- 
lar, the amplitude of the £ = I harmonic to the Fourier 
transform of the response function is proportional to a 
ppE theory parameter, which we are assuming here to be 
small, relative to the GR amplitude of the £ = 2 mode. 
Since GW detectors are much more sensitive to the GW 
phase than its amplitude, it is not entirely clear whether 
the £ =: I harmonic can be extracted for low signal-to- 
noisc ratio events. Whether the inclusion of such param- 
eter amplitude ppE parameters is of practical use can 
only be determined through an implementation of these 
proposals. 
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